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We examine the linear stability of static, spherically symmetric wormhole solutions of Einstein's 
field equations coupled to a massless ghost scalar field. These solutions are parametrized by the 
areal radius of their throat and the product of the masses at their asymptotically flat ends. We 
prove that all these solutions are unstable with respect to linear fluctuations and possess precisely 
one unstable, exponentially in time growing mode. The associated time scale is shown to be of the 
order of the wormhole throat divided by the speed of light. The nonlinear evolution is analyzed in 
a subsequent article. 
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> ■ I. INTRODUCTION 

o : 

Einstein's general theory of relativity leads to spectacular predictions like the big bang in the early universe and the 
Q\ . existence of black holes and gravitational radiation, all of which are subject to mainstream research in astrophysics 
t— I ' and cosmology. More speculative predictions of Einstein's theory are traversable wormhole geometries [l[ consisting 
of two asymptotically flat ends connected by a common throat through which time-like test particles might travel. 
• Similar constructions can be obtained by removing two holes from a given three-manifold and connecting one with 
Qi another by a handle, in which case wormholes have been proposed for interstellar travel or building time machines 

± : il- n 

j^Tjj However, as a consequence of the topological censorship theorem [J] wormholes must be supported by matter fields 
1—1 1 which are "exotic" in the sense that they violate the averaged null energy condition. In particular, this implies the 
' violation of the weak energy condition which means that some time-like observers measure a negative energy density. 
On the other hand, it has been shown recently [5] that it is possible to construct wormholes where the violation of 
the averaged null condition can be made arbitrarily small. Therefore, such wormholes could in principle be supported 
f^i , by quantum fields since quantum effects can lead to violations of the energy conditions. 

In this work we assume the existence of stationary wormhole geometries and analyze their stability with respect to 
small initial perturbations. If stable, the time evolution of such a perturbation may deform the wormhole throat by 
\^ • a small amount, but for large times the perturbation decays to zero and the equilibrium configuration is recovered. 
If on the other hand the wormhole is unstable, the time evolution may induce large deformations of the throat and 
in this case it is possible that the wormhole is destroyed making it unsuitable for interstellar travel or building time 
' machines unless the time scale associated to the instability is large enough. 

Clearly, the notion of stability depends on the precise matter model that is being considered. In this work, for 
simplicity, we focus on a massless ghost scalar field. By a ghost field we mean one which violates the dominant energy 
condition. A massless ghost scalar field consists of a minimally coupled, massless scalar field whose kinetic energy has 
a reversed sign. Such a field violates the null, and as a consequence also the weak and dominant energy conditions. 
We analyze the stability of static, spherically symmetric wormhole geometries supported by such fields. We prove in 
this paper that all such solutions are unstable with respect to linear perturbations and possess precisely one unstable 
mode which grows exponentially in time. Furthermore, we show that the associated time scale is of the order of the 
areal radius of the wormhole throat divided by the speed of light. In particular, this means that a wormhole whose 
throat has an areal radius smaller than lkm, say, decays in a few /is. 

Wormholes supported by scalar fields have been considered before (see, for instance, [f| 0, H, H, EH EH EH EH ) • ^ n 
[HSUI wormholes with a nonminimally coupled scalar field are shown to be unstable with respect to linear fluctuations. 
In Q wormholes with a massless, ghost scalar field are constructed. These wormholes are static, spherically symmetric 
and connect two asymptotic ends which either have both vanishing ADM masses or have ADM masses of opposite 
signs. Based on a linear stability analysis for the zero mass case it was claimed in Q that the wormholes with 
sufficiently small ADM masses are stable. However, a numerical simulation performed in indicates that small 
perturbations of the zero mass wormholes do not decay under nonlinear evolution and may cause the wormhole to 
explode or to collapse. Here, we point out that the linear stability analysis of 9] is incomplete in the sense that 
it only shows stability of the zero mass wormholes for a restricted class of perturbations where the areal radius is 
unperturbed and the perturbed scalar field vanishes on the throat. However, such perturbations are artificial in some 
sense because one could imagine perturbing the scalar field by a small ingoing pulse which is supported away from 
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the throat at some time t = 0, say. As t grows, this pulse travels towards the throat and since the wormhole metric 
is everywhere regular one expects it to reach and cross the throat at some finite time. On the other hand, requiring 
that the perturbed scalar field vanishes at the throat corresponds to placing a mirror at the throat which reflects 
the scalar pulse. As we show in this article, the zero mass wormholes are linearly unstable in the absence of such a 
mirror, i.e. if one allows perturbations which are not required to vanish at the throat. Furthermore, we generalize the 
stability analysis to the static, spherically symmetric wormholes with nonvanishing ADM masses and prove that they 
are linearly unstable as well. In a subsequent paper we analyze the nonlinear evolution of small initial perturbations 
by numerical means and show that depending on the details of the initial perturbation the wormhole either expands 
or collapses to a Schwarzschild black hole. 

This work is organized as follows. In section[TT]we review the static, spherically symmetric wormhole configurations 
which are supported by a ghost scalar field. In section HTT1 we perform the linear stability analysis. We do so by first 
generalizing the master equation obtained in [9( to massive wormholes. This master equation turns out to be singular 
at the throat for both massive and massless wormholes. As a consequence, it requires the solution to decay sufficiently 
fast as the throat is approached. We then transform the singular master equation in a regular one and show that 
the new master equation possesses an unstable mode which grows exponentially in time. The associated time scale is 
computed in section IIVI by a numerical shooting method. Finally, conclusions are drawn in section [V] 



II. STATIC, SPHERICALLY SYMMETRIC WORMHOLES SUPPORTED BY A GHOST SCALAR 

FIELD 

We consider a gravitational field g which is coupled to a massless ghost scalar field $. The action is 

16^G / (-K+^^-V^V^gd^x, 

where R and V denote the Ricci scalar and the covariant derivate, respectively, associated to g, and k = —8ttG 
with G Newton's constant. Notice that k is negative which is the reason for calling the scalar field "ghost". The 
corresponding field equations are 

R»u = kV,$-V„$, (1) 

= VV^S, (2) 

where R^v is the Ricci tensor associated to g. 

For a spherically symmetric field, local coordinates t, x, ip can be chosen such that 

g = -e 2d dt 2 + e 2a dx 2 + e 2c (d$ 2 + sin 2 i? dip 2 ) , (3) 

where the functions d = d(t,x), a = a(t,x) and c = c(t,x), as well as the scalar field $ = &(t,x) depend only on 
the time coordinate t and the spatial coordinate x. We are interested in traversable wormhole geometries which 
consist of a throat connecting two asymptotically flat ends at x — > oo and x — » — oo, respectively. This means 
that the areal radius r = e c is strictly positive and proportional to |x| for large |a;| and that the two-manifold 
(M,g) = (R 2 , —e 2d dt 2 + e 2a dx 2 ) is regular and asymptotically flat at x — > ±oo. 

For the spherically symmetric ansatz ([3]) the Einstein equations ([1]) yield the Hamiltonian constraint TL := 
~e a - d {R tt + e 2d R/2) + n(e a - d § 2 + e d ~ a $ 2 )/2 = 0, the momentum constraint M := -{R xt - n<$> t <& x ) = and 
the evolution equations e d ~ a (R xx — k^ 2 ) + H = and R$$ = R w = 0. This and the wave equation © yields the 



evolution equations 

d t {e a - d a t ) - d x (e d ~ a d x ) - e a ~ d c 2 + e d " a c 2 x - e a+d - 2c = ~ [e a " d $ 2 - e d ~ a <S> 2 x ] , (4) 

d t (e a ~ d+2c c t ) - d x (e d ' a+2c c x ) = ~e a+d , (5) 

d t (e a ~ d+2c <S> t ) - d x (e d - a+2c <S> x ) = 0, (6) 

and the constraints are 

H := e d ~ a [2c xx + (3c x - 2a x )c x ] - e a - d c t (2a t + c t ) - e a+d ~ 2c + | [e^ 2 + e d - Q $ 2 ] = 0, (7) 

M := 2c tx + 2c t c x - 2d x c t - 2a t c x + k = 0. (8) 



Here, the subscript t and x refer to the derivatives with respect to t and x, respectively. 
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For a static configuration, the scalar field $ and the metric coefficients d, a and c are independent of t. In this case 
the field equations can be integrated analytically [3, EH- In order to see this we adopt a gauge where d = — a in 
which case the field equations reduce to 



= 2, 



e 2 ^c, 



= 1, 



3 2(c-a) $a; 



= 0, 



(9) 



and 



2a x c x -e 2 ^-U 2 x = Q. 



(10) 



The first equation has the general solution e 2 ( c ~ a > = x 2 + 2ol\x + ao with two constants ao and a±. By a suitable 
translation of the coordinate x it is always possible to obtain ot\ — 0. Furthermore, since we are interested in wormhole 
geometries with the properties described below equation ([3]) we need e 2 ^'^ > for all x G E. Therefore, we have 
g2(c-o) _ x 2 _|_ ^2 w ith S0lne strictly positive constant b > 0. The second and third equations in ((9]) have the solutions 



c = — log(x 2 + 6 2 ) — 7i arctan ^— ^ 
with constants 70, 71, $0, It follows that 



7o , 



$ = $1 arctan ^-J + <t> , 



d = —a = 71 arctan(aj/6) + 70 



(11) 



(12) 



and the metric reads 



3 2 71 arctan(x/b)+2 70 ^2 + g -2 71 arctan^/f.) -2 7 o ^2 + ^2 + & 2) ^2 + ^2 ^ ^2)] 



(13) 



Notice that the constant rescaling t <— ► exp(— f2)i, x i— > exp(f2)a;, 6 i— ► exp(f2)6, 70 h ~ > 7o + ^> 7i ~~ * 7i with f2 
a nonvanishing constant leaves the metric unchanged. So for example, we can rescale the coordinates such that 
71 7r + 270 = in which case e 2a — > 1 and e 2c /x 2 — > 1 as x — > +00. Alternatively, we can rescale the coordinates such 
that — 71 7r + 270 = in which case e 2a — > 1 and e 2c /x 2 — ► 1 as a; — > —00. With these observations in mind we see that 
the metric ()13[) has indeed two asymptotically flat ends at x — > +00 and a; — ► — 00, respectively. Finally, equation (|10p 
yields the relation 



^ 2 = 2(l + 7l 2 ) 



(14) 



between the parameters 71 and $1, where we recall that k = —&t:G is negative. 

Summarizing, we obtain the bi-parametric family of wormhole solutions described by equations (|13p and (|11I12|) . 
This family of solutions has been obtained in [lj [TEI, and in Q in the context of wormholes. The solutions can 
be parametrized by the constants B := 6e~ 70 > and 71 > which are invariant with respect to the rescaling 
b 1— > exp(fi)b, 70 1—* 70 + discussed above. $1 is determined (up to its sign) by the constraint (fl"4|) . (The sign of 
$ is not important since the field equations are invariant with respect to the transformation $ 1— > — $. Similarly, 
the parameter <!>o has no physical meaning since only the gradient of $ appears in the equations. Furthermore, it is 
sufficient to consider nonnegative values for 71 since the solution is invariant with respect to a change of sign of both 
71 and x.) In order to relate the parameters B and 71 to physical quantities we first remark that the areal radius, 



has a global minimum at x = aJt/troerf = 7i& since c x — (x — ^ib)/{x 2 + b 2 ) and lim (r/|x|) = exp(— 70 =p 7 1 7r/2). 

X — >±0O 

Next, we compute the Misner-Sharp mass function [l6j |. For the spherically symmetric spacetime metric given by 
equation §5§ it is defined by 



T 6 

m (^ s ) == 2 [1 ~ -9( dr ' dr )\ = "9" 



1 + e 2( C ~d) c 2 _ e 2( C -«) c 2 



For the static solutions described by equations (|11I12[) this yields 

(x-716) 2 



m{x) = - 



1 



+ b 2 
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The ADM masses of the two asymptotically flat ends can be computed by considering the asymptotic values m±oo := 
lim m(x) of the mass function, which yields = -B71 exp(— 7!7r/2) and to_oo = — -B71 exp(7!7r/2). Notice also 

x — * : -x 

that at the throat, m(x t hroat) = r t hroat/^- 

Therefore, the parameters B and 71 determine the ADM masses of the two asymptotic ends, moo and m-00, and 

71 arctan(7i 



the areal radius of the throat, r throat = By/l + 7f 



Conversely, the areal radius of the throat and the 



product of the asymptotic masses uniquely determine the parameters B and 71. The particular case ji = yields a 
zero mass wormhole, m^c = m—00 = 0, in all other cases the asymptotic masses are nonzero and have opposite signs. 

Before concluding this section, we compute the violation of the averaged null energy condition according to [f|, 
which is given by the volume integrals of p + p? over the two asymptotically fiat ends, 



±00 



I± := 47r 



(p+Pr)r 2 r x dx, 



where p and pf are, respectively, the energy density and the radial pressure as measured by static observers. Since 
p = Pf for a scalar field and since the Hamiltonian constraint implies that 



(2m) x 



„3c-2a. 



= 



K\pr r x 



for static solutions we obtain 



I± = R (2- ± oo 



_ 8tt 

^throat) — "j 7 T throat 



±2- 



Ti 



v/1 + 7? 



7i[arctan(7i)=F7r/2] 



(15) 



We have /_</+< 87t|k| 1 r t ^ roat (2e 1 — 1) < 0, so both integrals are strictly negative as expected from the violation 
of the averaged null energy condition. However, we also see that both integrals can be made arbitrarily small by fixing 
71 and letting r throat — > 0. 



III. LINEAR STABILITY ANALYSIS 



In this section we analyze the linear stability of the two-parameter family of static wormhole solutions discussed in 
the previous section. In order to do so, we consider small perturbations of the form 

$(A) = $ + A<5$ + e>(A 2 ), 

where $ is the background solution, and where 

6<i> := A$( A ) 

denotes the variation of The same applies to the other fields d, a and c. Since there are no dynamical degrees of 
freedom for a spherical symmetric gravitational field we expect to be able to describe linear fluctuations by a single 
master equation for the linearized scalar field 5$. In fact this is known [ItJ to be the case for a large class of matter 
fields when considering linear fluctuations of static configurations. 

However, as we will see in Sec. IIII Al the derivation of the master equation for 5$ requires a gauge where the areal 
radius is unperturbed, i.e. where Sr — 0. Since r — e c is a scalar field on the two-manifold M its linearization Sr 
transforms as 



Sr h-> Sr + £ ? r = 6r + £ x r x (16) 

with respect to infinitesimal coordinate transformations on M which are parametrized by the vector field £. If r x =/= 
it is always possible to achieve the gauge Sr = by an appropriate choice of £ x . However, for wormhole topologies, 
r x = at the throat; hence it is not possible to set Sr — everywhere unless we restrict ourselves to perturbations 
which hold the areal radius of the throat (and any other surface of critical r) fixed. 

An alternative approach which does not require the gauge Sr = is to replace 8<$> with the linear combination 

* := e c \ 5$ - ^Sc^ , (17) 
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which is invariant with respect to infinitesimal coordinate transformations on M . In Sec. IIII Al we show that the 
linear perturbations can in fact be described by a master equation which has the form of a wave equation for with 
potential V. However, we will see that this potential diverges asi-> x t hroat approaches the throat. This property is 
not entirely surprising since for 5c ^ the gauge-invariant quantity \& is not well-defined at the throat where <& x /c x 
diverges. 

In Sec. IIII Bl we transform the singular master equation for ^ into a new master equation which also has the 
form of a wave equation with a new potential W which is everywhere regular, thus facilitating the stability analysis. 
Furthermore, we prove in Sec. IrnCl that the new equation possesses a unique unstable mode which grows exponentially 
in time. This mode is then shown in IIII Dl to give rise to an exponentially growing solution of the linearized field 
equations which is everywhere regular, thus proving the linear instability of the static, spherically symmetric wormholes 
supported by a massless ghost scalar field. 



A. Master equation for 9 



Here, we derive the master equation for \& which is valid for points away from the throat. In order to do so, we 
first consider the linearized Hamiltonian and momentum constraints which read 

8{e a - d H) = 2Sc xx + (6c, - 2a x )Sc x - 2c x Sa x - 2(Sa - 5c)e 2 ^ + K$ x 5$ x = 0, (18) 
5M = 26c tx + 2(c x -d x )6c t -2c x 6a t + K $ x 5<S> t = 0. (19) 

With the help of the background equations (e d ~ a+2c d x ) x = 0, (e d ~ a+2c c x ) x = e a+d and (e d ~ a+2c $ x ) x = 0, one can 
show that these equations possess the first integral 

Sc x + (c x - d x )Sc - cja = ~$ x 6$ + ae a - d ~ 2c , (20) 

where a is a constant whose meaning is explained below. 
Next, the linearization of equations (|5l6p yields 



Sett ~ e d - a - 2c d x (e d - a+2c Sc x ) = c x e 2(d - a) d x {5d - 5a + 25c) - 2e 2id ' c \Sa - 5c), (21) 
S$ u ~ e d - a - 2c d x (e d - a+2c 5$ x ) = <S> x e 2[d - a) d x {5d - 5a + 25c). (22) 

Now we are ready to obtain the master equation for For simplicity, we choose a gauge in which 5c = that is 
allowed since we are only considering points away from the throat. In this gauge, \& = e c <5$ and equation (|20|) reduces 
to 2c x 5a = K<S> x 5<&-2ae a - d - 2c . Equation jUJ then yields c 2 x d x (5d~6a) = 2c x c 2(a -^5a = ne 2 ^ a - c ^ x 5^>-2ae 3a - d - 4c 
allowing us to re-express the right-hand side of equation (|22[) in terms of (5$. The result is the master equation 



with the potential 



and the forcing term 



*« - e d ~ a d x (e d - a V x ) + V(x)« = Q(x), (23) 



vu) = <*■■'-> iv) ■ (24) 



Q(x) = -2a^e d+a - 3c . (25) 

C 



Since ^ has a gauge-invariant meaning, this equation holds for arbitrary gauges, and not only gauges for which 5c = 0. 
For the static wormhole solution (|llll2p we have 



gd— a arc t an ( a; /^)+27o 

„47i arctan(a;/6)+47o 

V{x) = - 



X 



+ b 2 



(.x- 7l 6) 2 ( 2(l + 7l 2 )& 2 ^ 



1 



(x - jib) 2 



g37i arctan(x/6) — 370 

Q(x) = -2cr$i6- 



Vx 2 + b 2 (x - 71 b) 5 
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We see explicitly that the potential V diverges at the throat x = -fib. Therefore, the master equation (|23[) requires 
^ to decay sufficiently fast to zero as x — > jib in order to be meaningful. If 6c = everywhere this in turn 
requires to be zero at the throat. As an example, the potential for the zero mass wormhole simplifies to V(x) = 
(b/x) 2 (3x 2 + 2b 2 ) I (x 2 +b 2 ) 2 and is positive away from the throat 1 . However, the two asymptotic regions are separated 
from each other by the infinite potential wall at x — which acts like a two-sided mirror. Therefore, it follows that 
with this mirror in place perturbations described by the field $ cannot grow in time. On the other hand, there is no 
physical reason for restricting ourselves to perturbations with zero $ at the throat. In the next subsection we show 
that equation (|23p can be transformed into an equation that is everywhere regular and that does not require \I/ to be 
zero at the throat. 

Before we proceed, we analyze the interpretation of the integration constant a. For this we notice that the lineariza- 
tion of the static expressions pi!12[) with respect to the parameters b, 70 and 71 must satisfy the master equation 
(|23[) for some value of a. Computing the corresponding variation, we find 

* = ( e*>8B - b + 7l ? X arctan(x/6)(57i^ , (26) 

x - "fib \ 1 + 7f / 

where B = be~ 10 . By inserting this expression into equation (|23p we find that 

a = B5-n + j r 5B = S(Bji). (27) 

Therefore, the integration constant a describes static perturbations which change the value of the product of the 
asymptotic masses, moom-oo = — (S71) 2 . Since any perturbation can be written as the sum of such a static pertur- 
bation plus a perturbation with S(B"/i) = 0, we can assume that a = for the following. The master equation J23|) 
then still admits the static solution 

*o = (l + x^T 2 ^ arctan^/fo)) , (28) 

x - 716 V b 1 +7i / 

which is obtained from (|26p with S^iB) = and 5B — — e~ 7o /<l>i. The existence of this solution plays an important 
role in the next subsection. 



B. Transformation to a regular equation 

Here, we transform the singular master equation ([23l into a regular one. For this, and for notational simplicity, we 
first rescale the coordinates t and x such that 6 = 1 and 70 = 0. Next, we define the two differential operators 

where d := e 271 arctan ( x ) Q x and Vl'o is the particular solution given in equation (|28| . Since A^A = —d 2 + (9 2 4'o)/4'o = 
— d 2 + V^q we can rewrite the master equation (|23|) as 

(d 2 + A^A) * - 0. (29) 

Applying the operator A on both sides of this equation, we find that the quantity x '■— A9 satisfies 

(d 2 + AA^) X = 0, 

where 

AA^-0 2 + d ^l = -8 2 -V + 2( d -^)\ 

Therefore, x satisfies the transformed equation 

Xtt ~ d 2 X + W(x) X = (30) 



1 The positivity of V led to the claim in [j| that there are no unstable, radial modes. 
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with the transformed potential W — — V + 2 f J . Explicitly, we have 

•y/J + a; 2 e _71 arctan ( ! >0 



*o(x) 



F(£c), := 1 + 



7i 



r(l + "fix) arctan(a;). 



x - 71 1 + 7i 

It is not difficult to prove that the function F is strictly positive. A short computation gives 



<9>F 



0_ p^7i arctan(a;) 



£k + x " 71 



1 



1 + X 2 IE — 7! 



Using this, a lengthy calculation yields the following expression for W: 



W(x) 



c 47i arctan(:c) 



3 I 3 (X " 71)2 I 2( F * 
l + x 2 (l + x 2 ) 2 \F 



47i F x 471 x - 71 



l + x 2 F 1 + 7 2 F(l + x 2 ) 2 



(31) 



(32) 



where 



71 



1 + T? 



1 ~T~ 'Yl^' 

— + 71 arctan(x) 

1 + x z 



F = 1 



71 



l + 7i 



(1 +7ix) arctan(x). 



For 71 = the new potential reduces to the simple expression W(x) = —3(1 + x 2 ) 2 which is strictly negative. For 
71 7^ we have 



W(7i) = e 471 arctan (Ti) 



27? 



1 + 7? (l+7i 2 ) 2 



< 



and W(x) = e ±2ni1 \2x~ 2 + 0(x~ 3 )] for x — > ±00. Therefore, is positive for large |a;| and negative near the 

throat x = 71. In the next subsection we prove that the new master equation (|30[) possesses a unique unstable mode 
of the form 

X (t,x) = e pt xo(x), 

where xo is the ground state of the Schrodinger operator H := AA* = —d 2 + W with negative energy Eq — —0 1 . 
Notice that Xo{x) decays exponentially to zero as |x| — > 00. 

C. Existence of a bound state with negative energy 

Here, we establish the existence of a bound state with negative energy of the Schrodinger operator H := AA'. For 
this, we first introduce the new coordinate 



p{x) 



-2-yi arctan(y) 



dy 



which monotonically increases from —00 to +00 and satisfies d — e 2 7i arctan(x)^ _ q anc j jj m pjx = exp(+7i7r) 



x — >±oo 



Since W is real and bounded, H = —d 2 + W(p) defines a self-adjoint operator on the Hilbert space L 2 (R, dp) 
with domain consisting of the dense subspace D(H) :— H 2 (M.,dp) of functions whose zeroth, first and second order 
derivatives are square- integrable on R. Let us assume first that 71 7^ 0. In this case the function 



x - 71 



7i arctan(a;) 



v 7 ! + x 2 F(a 



(33) 



is an eigenfunction of H with zero eigenvalue since W = ®od;(% L ) and since F grows linearly in \x\ for large |a;|. 
Additionally, this zero mode has precisely one zero, namely at the throat x = 71. It follows from the nodal theorem 2 



2 See, for instance Ref. flU . 
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that H possesses exactly one bound state \o with negative energy Eq = —0 1 < 0. Its value will be computed in the 
next section via a numerical shooting method. 

For 71 = the function l/^o = x/yl + x 2 is not an eigenfunction of H anymore because it does not decay to zero as 
|x| — > oo. Nevertheless, we may still expect the existence of a unique bound state with negative energy by continuity. 
In order to obtain an upper bound for the ground state energy E for the case 71 = where W = —3/(1 + a; 2 ) 2 , we 
use the Rayleigh-Ritz 3 variational principle, 

E = inf (34) 



where (• , •) denotes the scalar product of i 2 (M, dp). Testing with the particular functions ^ k{x) = (1 
K > 1/2, gives 



-K/2 



*k) = 2/jc, (tf H^ K ) = 2 [K 2 I K+l - (K 2 + 3)I K+2 ] 



where the integral 



IK 



dx 1 
K > - 







{l + x 2 ) K ' 2 



satisfies the recursion relation 



Therefore, we obtain the bound 



lK+1 = ~2K~ lK ' K> 2 



Ea< _l liK) , m ^-W;"-*> , K> 1. 

The function / is negative for 1/2<X<3 + 2V3 and positive for K > 3 + 2\/3- Its minimum lies near K = 2 for 
which f(K) = —11/2. Since the minimum of W is —3 we obtain the estimate —3 < Eq < —11/8 = —1.375 for the 
ground state energy. This estimate is used as a starting point for a numerical shooting method described in Sec. IIVI 

D. Existence of an exponentially growing mode 

To prove that the unstable mode found in the previous subsection gives rise to an exponentially growing regular 
solution of the linearized field equations we still need to show that there exist perturbation amplitudes Sd, 5a, 8c and 
5$ which satisfy the linearized field equations, are everywhere regular and grow exponentially in time. Here, we show 
that such a solution exists indeed. 

In order to do so we first notice that $ = Jv\ satisfies the master equation (|23l) . Next, we choose the infinitesimal 
coordinates such that 54> = 0. This is possible because <S> X = $i/(l + x 2 ) ^ for all x € R. In this gauge, 
$ = — &ie c 5c/(x — 71), and so we obtain 



eCSc = -^r 1 ^ = ^r 1 (d + ^ )x = - — 



2 



(.x - jx)8 x + (x- 71)— + 2 - 1 



(36) 



which is regular everywhere and decays exponentially fast to zero as |x| — > 00. We may obtain (5d and 5a from this 
using the perturbation equations (j2"D|) and (f2"2"| which, in the gauge 5$ = yield 

c x 5a = fcr + (ca - d x )Sc, 
5d x = 5a x — 2Sc x . 



3 See, for instance Ref. [ill . 
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A careful calculation using the fact that x satisfies (d 2 — W)\ — P 2 X an d the identity 

7i 1 7i0-7i) 



Fx 
F 



yields 



8a = 
Sd = 



Ti 



1 



2 1 



F(1 + 7i 2 )(1 + :e 2 ) 



[1 + x arctan(x)] 



x arctan(x) 



-1 



■7? 

9 

7i 



F 

x arctan(ir) 



Sc + 



vT 



7? 



6c + 



v / TT~ 



-71 arctan(a;) o2 



Fx, 



'-e-^ iaictan{x) (3 2 x + h(t), 



(37) 

(38) 
(39) 



where h(t) is an arbitrary constant which depends on t only and could be eliminated by an infinitesimal coordinate 
transformation. It can be checked that the expressions (|36|38|39|) also satisfy the evolution equations (|2"Tj) and the 
linearization of equation ^ with = 0. 

Therefore, we obtain an exponentially growing solution of the linearized field equations which is everywhere regular 
and decays exponentially to zero for \x\ — ► 00. This proves the linear instability of static, spherically symmetric 
wormholes supported by a massless ghost scalar field. Their nonlinear evolution are studied in a subsequent paper. 



IV. TIME SCALE OF THE INSTABILITY 

In this section we compute the ground state energy E Q = —f3 2 of the Schrodinger operator H via a numerical 
shooting method. For this, we first multiply the differential equation H\ = —fi 2 X 011 both sides by the factor 

e -4 7l arcta„( 7l ) and obtam 



-{a(x)d x } 2 X +[P 2 + W(x)]x = 0, 

where a{x) = e 2 7l[arctan(>)-arctan(7i)] ^ _ g-27i arctan(7i ) ^ an( j 

2 



(40) 



W(x) = a(xY 



1 



, (x ~ 7i) 
'(l + z 2 ) 2 



F x 



471 



x - 71 



4 7 i F x 
1 + x 2 F ' 1 + 7 2 F(l + x 2 ) 2 



For 71 — the potential is strictly negative. For 71 > it is convenient to rewrite the potential as 



W{x) = a{xf 



( 7l 2 -3-6 7 ix)/(x) 2 -27r 
(l + x 2 ) 2 l(x) 2 



1 7 2 (x- 7l ) 
F(l + 7l 2 )(l + x 2 ) 



l(x) 



7i 



(1 + x 2 )l(x) 



where l{x) := 1 + xarctan(a;) > 1 and where we have used the identity (|37[) . Since the second term inside the curly 
brackets is a perfect square we see from this representation that W{x) is positive for all values of x small enough such 
that 



(6 7 ia; + 3-7 2 )Z(a ; ) 2 < -2 7 f 



(41) 



Therefore, for small x, p 2 + W(x) > which means that x cannot have a local maximum (minimum) if x > (x < 0). 
For large x, on the other hand, W(x) — e 2 7i(T-2arctan(7i)) [2x~ 2 + 0(x~ 3 )] decays like 1/x 2 and so the normalizable 
solutions behave like 



X(x) 



,-/3p(a 



(42) 



for large x, where p(x) = f* 

Our shooting procedure consists in integrating the differential equation (|40p starting with the asymptotic solution 
(|42[) at some large value of x and integrating numerically towards small values of x. The integration is stopped as 
soon as the inequality (|4*T|) is satisfied and x ' Xx < 0. If this is the case, the solution must diverge as x — > —00 since 
\x\ cannot have a local maxima in the region where the inequality (|4ip is satisfied. Therefore, we look for the value 
of P for which x -> as 1 -> —00. 

For the numerical integration we compactify the domain by means of the transformation x = tan(7rz/2) which 
maps z € (—1,1) onto x S K and use a fourth order Runge-Kutta integrator with fixed step size. We start with 
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the massless case 71 = where the estimates in the previous section provide the bound a/11/8 < (3 < a/3 for the 
parameter 0. For /3 = a/11/8 the final value for \ aTi small x is negative, for /3 = a/3 it is positive. The optimal 
value for f3 is obtained by bisecting the interval [a/11/8, a/3]. Then, we increment the value of 71 by a small amount 
and repeat the above shooting method, finding the optimal value for /3 for the new value of 71. Continuing this 
way we obtain the results summarized in table [I] and figure [TJ where we show the associated time scale in terms of 
proper time at the throat of the background solution, r = e 71 arctan (Ti)^ ) anc [ units of the areal radius of the throat, 
rthroat = &a/1 + 7 2 e~ 71 arctan (7i). This time scale is given by 

271 arctan(7i ) 

c ' throat ' throat 

Tunstable — , „ n — ,-> ~= ■ 

a/1 + 7? P a/1 + 7? P 

As is apparent from the results in table [T] and figure [TJ the time scale converges to a value of 0.590r t ^ roat for large 
values of 71. The limit 71 — > 00 is analyzed next. 



71 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.2 


1.4 


1.6 


1.8 


2.0 


T 


0.846 


0.841 


0.826 


0.805 


0.782 


0.758 


0.737 


0.718 


0.701 


0.687 


0.675 


0.656 


0.642 


0.631 


0.624 


0.618 


7i 


2.5 


3.0 


3.5 


4.0 


4.5 


5.0 


5.5 


6.0 


6.5 


7.0 


7.5 


8.0 


8.5 


9.0 


9.5 


10.0 


T 


0.608 


0.602 


0.599 


0.597 


0.595 


0.594 


0.593 


0.592 


0.592 


0.591 


0.591 


0.590 


0.590 


0.590 


0.590 


0.590 



TABLE I: Numerical values for the dimensionless time scale T := T unata bie/ rthroat for different values of 71. 




FIG. 1: The dimensionless time scale T := T unsta bie/ rthroat vs. 71. This plot and the values in the previous table strongly 
suggest that T converges to the value of 0.590 as 71 — > 00. The dashed line represents the value for T at 71 — > 00 which is 
computed in the next subsection. 



A. The large 71 limit 

In order to understand the behavior of the time scale for large 71 we first notice that the inequality (j4T|) is satisfied 
for x < 71/6 — 1/(271) — l/(37i) =: ^(7x)- Therefore, the potential W(x) is strictly positive for all x £ (—00, i?(7i)) 
and the solution grows at least exponentially fast with increasing x on this interval. Since i?(7i) — > 00 as 71 — > 00, 
the solution must vanish for each fixed x as 71 — > 00 if it is required to have a fixed, finite value at x — R{^\). For 
this reason, we first replace the coordinate x by z :— x/71 in the following, and then take the limit 71 — > 00. 

In terms of the new coordinate z the differential equation (j40|) is transformed into — [a(z)d z ] 2 x + [0 1 + W]x = 0> 
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Where Cl(z) = e 27l[arctan(7i2)-arctan( 7 i)] ^ ^ = ^ &nd 



W(z) = a{zf 



3" 



, , 7?(z-l) 2 

l + 7i 2 ^ 2 (l+7i 2 ^ 2 ) 2 



47? F, 



H - 1) 



l + 7l 2 z 2 F l + 7 1 2 F(l + 7l 2 z 2 ) 2 



where F z = 7 2 (1 + 7 2 ) 1 [(1 + 7 2 z)/(l + 7 2 z 2 ) + 71 arctan(7iz)] and F = 1 + 71(1 + 7 2 ) 1 (1 + jfz) arctan(7iz). From 
the above, W(z) > if z < 1/6 - l/(27 2 ) - l/(3 7 f). Now let z > be fixed. In the limit 71 -> 00 we have 



F 
7i 



71 



7T 

2 



and 



71 [arctan(7iz) — arctan(7i)] = 71 arctan 
Therefore, the differential equation (|4"0)) converges pointwise to 
- (V^"^) 2 X + e 4(1 - z ~ 1} [2z- 2 - 10* 



-3 



7i(*~l) 
1 + 7l 2 z 



3aT 4 ] X 



z- 1 



-/3 2 X, 



z > 0. 



(43) 



The operator on the left-hand side is a formally self-adjoint operator on the Hilbert space L 2 ((0, 00), e _2 ^ 1_z ^dz) 
which is singular at z = 0. We obtain a zero mode of this operator by multiplying equation (|33ll by 7ie~ 71 arctan (7i) 
and taking the limit 71 — > 00. This yields 



Tl ^ — 71 arctan(7i ) 7l jv 7l [arctan(7iz) — arctan(7i)] 

*0 ~ y/l + 7l 2 ^ 2 ^ 



2 z - 1 



,1-2" 



(44) 



Indeed, it can be checked that fH)) solves the differential equation (jlB")) with /3 = 0. Furthermore, this solution decays 
like 1/z for large z and as z — » it decays rapidly to zero. Since it has exactly one zero, it follows again from the 
nodal theorem the existence of a unique bound state with negative energy Eq = —p 1 - A numerical shooting procedure 
similar to the one described above which starts at some large value of z, where x(z) w e ~Pp( z ) ^ p(z) = e -2 ' 1 ~'» ^dC, 
yields 1//3 = 0.588. Since on the other hand T unsta bie / nhroat — l/(yi + 7 2 /?) = 1/(3 in the limit 7! — > 00 this matches 
well the asymptotic value obtained in table [U 



V. CONCLUSIONS 



In this article we analyzed the stability of static, spherically symmetric general relativistic wormhole solutions 
sourced by a massless ghost scalar field. We found that all these solutions are unstable with respect to linear 
fluctuations of the metric and the scalar field. Furthermore we showed that the time scale associated to this instability 
is of the order of the areal radius of the throat divided by the speed of light, which is of the order of a few microseconds 
for a throat of radius of the order of kilometers. Therefore, the instability we found is likely to introduce a rapid 
growth or collapse of the throat which eventually may destroy the wormhole. In a subsequent paper [lj| we follow 
the nonlinear evolution of the unstable mode by numerical means and show that the wormhole either expands rapidly 
or collapses to a Schwarzschild black hole. 

Our conclusion about the instability of the zero mass wormhole is different from the results presented in Ref. 0]. 
As explained in Sec. IIII Al this is due to the fact that the results in Q only apply to a restricted class of perturbations 
which are required to vanish at the throat. As shown in this article, both the massless and massive wormholes are 
linearly unstable for the more general class of spherically symmetric perturbations which do not necessarily vanish at 
the throat. 

The question arises whether or not there exist wormhole solutions different from the ones considered in this article 
which are linearly stable. Among the potential options are: i) wormholes with a single asymptotic end obtained by 
removing two holes from M 3 and connecting them to each other by a handle, ii) rotating wormholes, iii) wormholes 
which are supported by a different kind of matter field, iv) wormholes in modified or higher-dimensional theories of 
gravity. Regarding option i), it is conceivable that a change in topology affects our instability result. Indeed, since 
the linear stability of a given solution is directly related to the spectral properties of the linear operator appearing 
in the perturbation equation it depends on the global structure of the wormhole, and not only on local properties of 
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the throat. Concerning option ii), the possibility that the rotation could stabilize a wormhole has been explored in 
where an exact solution describing the inner region of a rotating wormhole is constructed. However, based on our 
instability results, our expectation is that the centrifugal force due to rotation only shifts the equilibrium between 
the attractive gravitational force and the repulsive force of the ghost scalar field, but does not change its stability. 
At least we expected that by continuity, slowly rotating wormholes are unstable. Regarding options iii) and iv) there 
is already a vast amount of results in the literature. To mention only a few examples, it has been shown that static 
spherically symmetric wormholes in scalar-tensor theories of gravity are linearly unstable @, H| . On the other hand, 
a class of wormhole solutions which are supported by phantom e nerg y have recently been constructed in (20l . l2l| and 
shown to be stable with respect to perturbations inside this class [22|. 
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